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Abstract 

Within the framework of the Dyson-Schwinger equations in the axial gauge, we study 
the effect that non-perturbative glue has on the quark propagator. We show that Ward- 
Takahashi identities, combined with the requirement of matching perturbative QCD at 
high momentum transfer, guarantee the multiplicative renormalisability of the answer. 
Technically, the matching with perturbation theory is accomplished by the introduction of 
a transverse part to the quark-gluon vertex. We show that this transverse vertex is crucial 
for chiral symmetry breaking, and that massless solutions exist below a critical value of 
as- Using the gluon propagator that we previously obtained, we obtain small corrections 
to the quark propagator, which keeps a pole at the origin in the chiral phase. 
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1. Introduction 



The Dyson- Schwinger (DS) equations of motion are one of the main tools for the investiga- 
tion of non-perturbative effects. These equations are particularly suited to the study of the 
interface between perturbative and nonperturbative regimes, as they allow us to compute 
the evolution of Green functions when one enters the infrared region, and hence suggest 
modifications to the usual Feynman rules. 

In a previous paper, we have studied |I| the DS equations in the quenched approxima- 
tion, i.e. for pure gauge QCD, and found that, as was first suggested by Cornwall 0, some 
solutions for the gluon propagator are flatter than a pole in the infrared region. In order 
to avoid the question of ghost propagators, we worked in the axial gauge. Other groups 
H, §, 1^ have also found that flatter solutions exist, in covariant gauges, although it is at 
present difficult to relate our results to theirs. The absence of a pole in the gluon propaga- 
tor at fc^ = is natural if one assumes that gluons do not propagate to infinity, i. e. these 
solutions should describe confined gluons. Furthermore, as was pointed out by Landshoff 
and Nachtmann , the existence of such solutions is highly desirable in phenomenological 
applications, and gives us a practical way to extend usual perturbative estimates to the 
strongly-interacting sector of the theory. The use of these solutions has already met with 
some success in diffractive calculations 0. 

The DS equations constitute an infinite tower of integral relations between n-point 
functions. Hence, by themselves, they cannot be solved. For instance, if we consider the 
DS equation for the quark two-point function (the propagator), it involves a three-point 
function (the quark-gluon vertex) as well as the gluon two-point function (the propagator). 
One can imagine solving the pure-gauge DS equations for the gluon propagator as a first 
approximation, hence we shall assume here that the gluon propagator is known. However, 
resorting to higher-order DS equations to calculate the vertex will only bring in higher- 
order n-point functions, which are also a priori unknown. Hence one needs to use another 
input. 

The simplest assumption is to take the perturbative 7^ vertex, which is the so-called 
"rainbow" or "ladder" approximation. One can be however more sophisticated, and follow 
the observation of Baker, Ball and Zachariasen (BBZ) [|], namely that one needs to choose 
a vertex that will obey the Ward- Takahashi-Slavnov- Taylor identities. These constitute 
another nonperturbative statement of field theory, and hence must be valid in general. 

In our study of the gluon propagator [|1| , we chose for the three-gluon vertex the simplest 
function that would obey the Ward-Slavnov- Taylor identities. We showed that besides the 
original 1/k^ solution, there exists another solution, behaving like l/k^'"^ near k"^ = 0. To 
obtain a full picture of nonperturbative effects in QCD, one then needs to calculate the 
quark propagator corresponding to that gluon propagator. As we shall see, the problem of 
multiple solutions arises again, as for sufficiently small 0:5, a chiral solution exists together 
with a massive one. We shall show that one can constrain the spectral density of the gluon 
propagator and the value of the coupling so that only the massive solution survives. 



As the equation for the quarks is much simpler than that for the gluons, we investigate 
in more detail the effect of the ansatz for the quark-gluon vertex on the quark propagator. 
Namely, the minimal vertex, used by BBZ for gluons and by Ball and Zachariasen (BZ) for 
quarks |]T2[|, leads to the breakdown of multiplicative renormalisability. We demonstrate 
that, following the method of Curtis and Pennington (CP) 0, it is possible to recover 
it by choosing a specific form for the quark-gluon vertex, which agrees with perturbation 
theory in the ultraviolet region, and which generalizes the B(B)Z ansatz. We also show 
that massive solutions exist only for the CP vertex. Simpler ansatze lead to an inconsistent 
ultraviolet behaviour, unless the propagator is massless. We also explain under which 
conditions massless solutions will exist, even for the CP vertex. We proceed to solve the 
DS equation in the massless phase, and show that nonperturbative effects do not remove 
the pole of the quark propagator in the chiral phase of the QCD vacuum, suggesting that 
confined quarks can exist only after chiral symmetry breaking. 

This paper is set out as follows. In Section 2, we review our results for the gluon 
propagator, and outline the formalism that lead to these. We point out the problems linked 
with the BBZ choice of a three-gluon vertex. In Section 3, we discuss the DS equation for 
the quark propagator in the axial gauge, and the approximations made when imposing the 
Ward-Takahashi identity. In Section 4, we consider the part of the propagator that preserves 
chiral symmetry. We show that the problem can be reduced to a one- dimensional integral 
equation, which we then renormalise. We show explicitly that we recover multiplicative 
renormalisability. In Section 5, we give a criterion for the existence of massless solutions, 
and then proceed to the numerical solution of the equation in the massless case, for the 
gluon propagator that we previously derived. 



2. The Gluon Propagator 

The fundamental observation [^] is that one can use the Ward-Slavnov- Taylor (WST) iden- 
tities to obtain a closed equation from the DS equation, if one neglects quark loops and 
works with gluons only. In the axial gauge {n.A = 0), where we do not have to worry 
about ghost degrees of freedom, the DS equation relates the propagator to the three- and 
four-point vertices. 

The three-point vertex, F^^^ can be split into a part P^ transverse to the external 

gluon momentum A;^, Pj. ■ k = 0, and a longitudinal part, P]^ , with P]^ ■ A; 7^ 0. As 
the vertex is a third rank tensor, depending on the three vectors entering the vertex, this 
spht is not uniquely defined. The WST identities relate P]^ to the propagator, and these 



identities can be "solved", hence producing an ansatz for the longitudinal part 0, |T2]. It 
is this ansatz that we define as being the longitudinal part of the vertex in the following. 

One can then obtain a closed equation for the propagator if one makes two assumptions. 
First, the full propagator is supposed to have the same spin and gauge structure as the free 
propagator: 



.2^, _ nO (1. ^\^\^ 



D,,{e) = -D%{k,n)^ (2.1) 

where 

V^,{k, n) - V —^ + (^.^)2 l^-^) 

This spin structure projects out the four-point vertex contribution to the equation. If one 
further assumes that the three-point vertex is dominated by its longitudinal part F^^ , which 
is itself known from the propagator via the WST identities, one obtains an equation which 
involves only the gluon propagator. 

This equation still needs to be renormalised. One can show without approximation § 
that a propagator with the spin structure ( |2.1| ) has to be singular as k"^ -^ 0. Thus the 
inverse propagator should vanish at the origin. This allows one to subtract the quadratic 
divergences, thus renormalising the gluon mass to zero and cancelling the tadpole graphs. 
One is then left with logarithmic divergences, which can be dealt with through wave func- 
tion renormalisation at a renormalisation point ^g, by defining 

Zik') = Zif^DZnik') (2.3) 

However, the resulting definition of the renormalised coupling constant ag is slightly dif- 
ferent from the usual one: 

"^^^^^ 1 + a,Z{f^l) [/r rfPvl/(P, Z)]^^^^ ^ ■ > 

where a^ is the bare coupling and ^ a function linear in Z, the explicit form of which 
can be found in ref. [|ll]. It has to be noted at this point that Eq. ( p. 41) leads to two 
related problems. First of all, the propagator asymptotically behaves like l/log(fc^)^^/^^ 
as fc^ ^ oo [H, which is not what is obtained in perturbation theory. Furthermore, one 
obviously loses multiplicative renormalisability. This has the nasty consequence that, if 
the same truncation is applied to the quark propagator, the renormalised couplings are not 
equal anymore, because the functions \E' are different for quarks and gluons. 

The solution to the resulting equation, found in ref. 0, has a very singular behaviour 
proportional to l//c^ as the momentum squared goes to zero. This was advertised as a signal 
for a confining potential. It has to be pointed out that in order to renormalise the equation 
in the presence of a 1/k^ singularity, one needs to perform additional subtractions, and 
hence one obtains such solutions only if one assumes they exist. In ref. [|l], we exploited 
the fact that the integral equation obtained in ref. [^ is non-linear and can in general 
admit several solutions. We solved the original equation, without additional subtractions, 
and obtained a gluon propagator which, although singular as fc^ ^ as required in the 
axial gauge, has a cut singularity, flatter than a pole. This describes a soft gluon which 
is confined rather than confining and can conveniently be used in models of the pomeron 
such as that suggested by Landshoff and Nachtmann B . This propagator may be written 



in the axial gauge {n.A = 0) as 



A;2 



f^n 



h.2 



+ «2(^ 



62 



+ c In 



d{^ 



(2.5) 



The dimensionless constants ai, 61, 02, &2, c, (i and e, are given in the Table. The mass scale 
/ig is the value of momentum at which ZR{k'^) takes the value 1, so that the propagator 
matches the value of the free gluon propagator. It is not obtained from the DS equation, 
which is scale invariant, but rather from a fit to the total cross-section for proton-proton 
scattering calculated using this propagator. An optimum value of 0.8 GeV was found for 
/ig. It is worth pointing out that most of the ideas discussed here can be applied to other 
theories, living at other mass scales, such as technicolour, and that one would only need to 
adjust II g to do so. 



Table: The constants used in Eq. (pT 



ai 


0.88 


c 


0.59 


02 


-0.95 


d 


2.1 


hi 


0.22 


e 


4.1 


b2 


0.86 


agifig) 


1.4 



Note that in the case of a multiplicatively renormalisable theory, the actual value of as 
should not matter. In principle, one could always rescale as by varying /i according to 



asilJ,')Zji{fi'^, q^) 



«s(/w)-2R(/i^ g^) 



(2.6) 



where the first argument of Z refers to the renormalisation scale: Z{fi'^,fi'^) = 1. However, 
our gluonic solution does not obey these equations because of ( p.4|) , and one obtains a 
solution only for a given valueQ This value of ag{fig) enters phenomenological calculations, 
which then determine fig. 



Hence the assumption of a given spin structure (|2.1|) and the neglect of the transverse 
part of the three gluon vertex leads to a DS equation which admits two solutions, and 
which breaks multiplicative renormalisability. In order to explain the discrepancy of the 
BBZ result with more recent ones, the assumed spin structure of the gluon propagator has 
been put into question ||^. However, there is no doubt that such solutions exist. Further 
solutions can of course exist, given that the equation is highly nonlinear. The neglect of 
the transverse part of the vertex, on the other hand, may have important consequences, 
which we shall explore in this paper for the quark propagator. 



'^Note that there was a mistake of a factor 2 in our definition of ag in Refs. m 



3. The Dyson-Schwinger equation for quarks 

Following their work with Baker on the gluon propagator ^, Ball and Zachariasen have 
considered the DS equation for the quark propagator [jl2|. Choosing again a purely longi- 
tudinal quark-gluon vertex, they found that the quark propagator corresponding to a 1/fc^ 
gluon was suppressed, and became constant near the origin. Furthermore, they showed 
that a chiral solution was always possible, and argued that a solution could be found 
which would break chiral symmetry. We shall see that, in the case of a less singular gluon 
propagator, only the chiral solution is possible if one neglects the transverse part of the 
vertex. However, implementing the improvements proposed by Curtis and Pennington 0, 
we shall show that multiplicative renormalisability (and hence the transverse part of the 
vertex) leads to an equation that allows chiral symmetry breaking, and even suggests a 
range of parameters (describing the coupling and the gluon propagator) which would lead 
to a unique, massive, solution. 



3.1. The Ball-Zachariasen equation 

We define the propagator for a quark with momentum q as: 

%) = Fiq'h.q + G{q') (3.1) 

so that F{q'^) represents the chiral-symmetry conserving part and G{q^) represents the 
chiral-symmetry breaking part. We shall also use the equivalent notation: 

The free propagator is therefore obtained by setting F(g^) to 1/g^ and G{q^) to zero 
(neglecting all current masses). In the axial gauge Eq. (|3.1|) is not the most general form for 
the quark propagator. First of all, the functions F and G can depend on n.q as well as q^. 
As in the case of the gluon propagator P, [l| , we seek solutions for which these functions are 
independent of n.q. Furthermore, the gauge dependence can also arise through extra spin 
structures, proportional to 7 ■ n and n.q, which although absent for the free propagator can 
in general occur for the dressed propagator. As pointed out in ref. [0]), these extra terms 
drop out of the equations for F and G if one specializes to a gauge vector orthogonal to q. 
This gauge choice n.q = is the only one that makes the algebra sufficiently tractable that 
reliable solutions to the DS equations can be obtained. The question of the dependence of 
the quark propagator on n, although important, is difficult to address in the axial gauge, 
and beyond the scope of this paper. 

The DS equation for S'(g), sketched in Fig. 1, in Euclidean space is given by 

1 =^.qS,{q) - CFakJ^J^Zk{e)rD,,{k,n)Sb{q-k)T''{q-k,q)Sb{q) (3.3) 



Figure 1: A pictorial representation of Eq. (^.3|). The hatched circles represent the exact two-point 
functions, and the cross-hatched circle the exact three-point function. 



where the subscript b on a, Z and 5* indicate that these are bare quantities which will have 
to be renormalised. Cp = 4/3 is the quark Casimir invariant. 

In general it is not possible to solve this equation because it involves the unknown quark- 
gluon vertex function r'^(g — k,q). This vertex function is related to a four-point Green 
function via a DS equation. Thus we get an infinite tower of coupled integral equations. 
On the other hand, in the axial gauge, this quark-gluon vertex function connecting a quark 
of momentum p = q — k to one of momentum q through the absorption of a gluon of 
momentum k, obeys the Ward-Takahashi (WT) identity sketched in Fig. 2: 



rr,(p,g) = S{q) 



-1 



S{p) 



-1 



(3.4) 



Figure 2: A representation of the Ward-Takahashi identity of Eq. (^^) using the same convention 
as for Fig. 1. 

One can then solve this equation to determine F^, the part of the vertex longitudinal 
to k. Following ref. [|T^, we obtain: 



-S{q)T^(p,q)S{p) 



+ 



mp) + F,{q)]Y 



- \F,{p) - F,{q)] ^,_^2 



mp) - F,{q)] 



q2 _^p2 



r 



[Gtip) - G,(g)] 



7^ 7 ■ p + 7 ■ g 7^ 

p2 _ g2 



(3.5) 



The simplest approach consists, as suggested in ref. ||T3[, in neglecting the transverse part 
of the vertex: one then assumes that the integral equation is dominated by the part ( |3.5| ) 
which is determined from the above WT identity.^ The resulting integral equations for 
the two functions Fi,{q^) and Gi,{q^) separate and we have 



1 = q^F,{q')-Cpa,J^;^^^Me)D,,{k,n 



' F,{q') + F,{p' 



-s. 



flU 



+ 



(g^ — p2 



Fbiq ) - Fb{p ) -k 6^^ + qf,pu + Pf,q^ 



(3.6) 



and 



= Gb{q')-CFabJ^J^Zb{k')D,,{k,n 
,2\ ry /„2 



(g^ — p'^) 



Gbiq ) - Gb{p ) (g ■ k5f,„ + q^p^ + p^q^) 



(3.7) 



^In other gauges the Ward identity Eq. (3.4) must be replaced by the full Slavnov- Taylor identities |1J] 
which involve the interactions of Faddeev-Popov ghosts. Nevertheless it was argued in ref. |lq| that the 



contributions from Faddeev-Popov ghosts are small and so the Ward identity of Eq. (3.4) was imposed on 
the integral equation formulated in a covariant gauge. 



where again we have introduced the subscript b on F and G to indicate that these quantities 
are to be renormahsed. (There is no renormahsation for E since this function is zero for 
the free propagator, i.e. in the absence of an exphcit current mass there is no parameter 
in the QCD Lagrangian to renormahse, hence the ratio F/G does not get renormahsed). 



3.2. Multiplicative renormalisability and the transverse part of 
the vertex 



As was observed in [|T2|, Eq. (|3.6| , |3.7| ) are not multiphcatively renormahsable. In other 
words, if we renormahse the quark and gluon wavefunctions by imposing Fb{q^) = Fb{ii'j) 
Fulq"^), Gb{q^) = Fblfij) G^^q^) and Zi,{q^) = Zi,{^'^)ZR{q'^), additional terms need to be 
introduced in the definition of the renormahsed couphng, as in Eq. (|2.4| ). These terms 
are not the same in the quark and in the gluon case, so one loses not only multiplicative 
renormalisability, but also the universality of the QCD coupling constant! 

These problems come from the ultraviolet region, and can be traced back to the neglect 
of the transverse part of the vertex in the solution of the WT identities: the vertex (|3.5| ) 
does not match the perturbative one at high momentum transfer. Hence, it is necessary 
to postulate a transverse part that will restore multiplicative renormalisability. Curtis 
and Pennington have shown [^ that this goal can be achieved in QED by considering the 
perturbative limit of the vertex. Their argument can be trivially extended to QCD in the 
axial gauge: the one-loop corrections to the propagator are identical to those of QED up to 
the quark Casimir invariant Cp, and the vertex corrections are also the same: because the 
vertex and the wavefunction renormahsation constants Zi = Z2 are equal in this gauge, 
and because Z2 is a function of Gp only, the diagrams involving the three-gluon vertex, 
which depend on Gp — Ga_/2, with Ga the adjoint Casimir invariant, have to cancel. Hence 
the argument is totally similar to that of ref. P] and goes as follows. 



In general, the quark propagator has the perturbative limit 

^^In-^ 

47r A' 



lim^(g2) = l + ^ln^ (3. 



with ^ the anomalous dimension. The one-loop vertex can be shown ^ to tend to 

as ^ 



hm TP;'\p,q) = 7m 



1- 




47r \A2 



(^r, In I -r^ I 

- -^ — V^(^'P^m7-P-9m7-P + '^9m7-p) (3.9) 

One can compare this answer with the vertex (|3.5| ) in the same limit, to conclude that 
the transverse part must give: 

,lip T {p,q) = — ^ {-q^,l■p + l■q1^,l■p) (3.10) 



This tensor is indeed transverse to g ~ A;. One can then extend the tensor structure of 
this vertex so that it becomes transverse to k for any value of p and q. We find that the 
simplest extension is: 

sipKip. ,)s(,) = (m - m) ^>tf-p')-fa>-+p>-)(^-^-^-rt (3.11) 

where V is an expression symmetric in p and q, and behaving as g^ in the large-g^ limit. 
T> must not introduce any singularity, must be symmetric in q and p, an must satisfy 
r^(p, p) = 0. In general, 

^-if^P^ij^^yJ^) (3.12) 

with 7] a regular function such that 77(0, 0) = 1. 



It is worth noting that the tensor structure of the transverse vertex ( |3.11| ) can be 



generalised. Indeed, Ball and Chiu DTOf have given a set of 8 independent tensor structures 



spanning the space of regular transverse vertices. Only three of these have the correct 
helicity structure to contribute to the DS equation in the chiral limit, and are (in the 
notation of Ball and Chiu): 

T2' = [p^{k.q)-q^{p.k)]^-{p + q) 

U = r{q'-v')-{q + vYl-{q-v) (3.13) 

The transverse vertex given by Eq. ( |3.11| ) corresponds to: 

T^{P, q) = l: f ^ - ^1 iW + ^')T. + W- q')T, + ■■] (3.14) 



where the ellipses refer to terms that vanish in the DS equation. This differs from the 
structure of ref. 0, where the transverse vertex was chosen to be proportional to Tg. This 
can be understood from the fact that the expression of the longitudinal part of the vertex 
used here differs from theirs, by a transverse tensor proportional to T3, in such a way that 
we get the same high-g^ or p^ leading term for the total vertex. The addition of such extra 
subleading terms is always possible, and we shall investigate the effect of these through the 
variation of function D in Eq. ( p.ll| ). 

The transverse vertex ( |3.11| ), together with the longitudinal one ( |3.5| ), can then be 
included in Eq. ( p. 3D to obtain: 

d^k 



p'F,{p) + CFa,j -^;^Dl,Zi,{k 



-2 t^''^^) + ^^(^)] ^^"' ''2,2' — \P^^^^ + ^M?;^ + 2 *^^ ~ q)^f^A 

{p^Fbjp) - q'^Fbjq)) ( V(g^ - p^) + p^,Pu - q^,qu + p^,qu - Puq^,) \ .^ ^^. 



4. Solving the equation 



4.1. The angular integral 



The angular integral in Eqs. ( p.6| , T!7, |3.15| ) cannot be performed analytically without 
prior knowledge of the function F, since this occurs with the argument (g — k^. The 
approximation, first proposed by Schoenmaker [0] in the gluon case, consists in replacing 
a function f{{q — fc)^) by /(g^ + /c^). This is clearly valid in the regions fc^ >> q^ and 
k'^ << q^. One then needs to find a function / which is sufficiently slowly varying that the 
error generated by the integration over the whole range of fc^ be negligible. 

The first choice would he f = F. In this case, the contribution of the longitudinal part 



of the quark-gluon vertex to the DS equation vanishes. Hence the BZ equations (|3.6| , [3?^ ) 
admit only the trivial solution F{q^) = 1/g^, G = 0. If the transverse part of the vertex 
is added, the resulting equation ( p.l5| ) leads to a propagator close to the perturbative one, 
with a pole at k"^ = 0. Neither solution is a slowly varying function, hence the approximation 
is not justified in this case. 

In the following, we shall use j = T . Our solutions show that this assumption is 
reasonable so that the above ansatz is justified a 'posteriori. The only region where this 
approximation may have led to substantial errors would be for small g^ ~ fc^, but then 
the integrand is suppressed because (as can be seen from Eq. (|2.5|)) the function Zji^k"^) 
vanishes. This approximation keeps the contributions both from the longitudinal and 
from the transverse part of the vertex, hence we shall obtain nontrivial solutions both 
for Eqs. (p], ^ and (|3J[5| ). 



With this approximation the angular part of the integration over momentum k in 
Eqs. ( |3l6| , pl7| , p.l5|) may be performed, and both equations can be recast into the fol- 



lowing form if G = 0: 
1 = q^F,{q^ 



dk^ 



Z{e)A^{e/q') F,{q') + Z{e)A,{k'/q')F,{k' + q^) (4.1) 



The kernels Ai and A2 are given, both in the case of Eqs. (|3.6|, 3.7) and in the case of 



Eq. ( ^.15[ ) in the Appendix. This is a Fredholm equation of the second kind which must 



now be solved numerically. Before we can do this, however, we must consider the question 
of the renormalisations required to absorb the ultraviolet divergences. 



4.2. Renormalisation 

The renormalisation procedure is similar whether one neglects the transverse part of the 
vertex or not. The renormalisation constant Zb{fi^) has been introduced in Eq. (p.3| ). 
As mentioned in Section 2 we choose this constant such that the renormalised function, 



10 



^ni.'f) takes the value 1 at g^ = /i^, with jj,g taken to be 0.8 GeV. We hkewise introduce a 
renormahsation constant J-'b{^f) such that 

and 



are ultraviolet finite and we choose it so that the quantity J-'rIq'^) takes the value 1 at 
q'^ = fJ''f- Eq. ([4.1|) then becomes finite if one rewrites it in terms of the coupling 

a, = "^^(^/) (42) 

^ l-abZ{i2f)CF/4:'Kjdk^Z{k^)Ai{k\l2}) ^'' 

and the renormalised equation ( [4.1[ ) becomes: 

1 = {l-^^Jdk'Zn{k')[A,{k\q')-A,{k\f,})]yR{q') 

- ^^jdeZn{e){A,{k\q')FR{e + q')-A,{k\^])FR{k' + ^^])) (4.3) 

1. Renormalised BZ equation 

In the case of the BZ kernels, the UV divergence of the equation is concentrated in the Ai 
term, which has a log singularity as /c^ — >■ cx3. This means that the coupling ( [4.2|) cannot 
be related to the usual renormalised QCD coupling 

asi^if) = Zb{nf)ab (4.4) 

Hence one loses the universality of the QCD coupling. 



Furthermore, the equation for the chirality breaking term G is identical to Eq. ( |4.3| ) 
with the replacement F ^ G, 1^0. The leading behaviour at high g^ comes from the 
remnant of the UV term. If we assume, in agreement with the renormahsation group, that 
Zniq"^) ~ l/[g^log(g^)] as g^ — >• oo, the leading terms of the equation for the G term are: 

^ Gniq') - ^^ log(log(g2))G^(g2) (4.5) 

47r 

which has only G = as a consistent solution. Hence, unless the gluon propagator is highly 
singular, the BZ equation does not lead to chiral symmetry breaking. 

If one overlooks these problems, and goes ahead to solve the equation, one obtains a 
function -FR(g^) which is very close to the perturbative 1/g^. The propagator in this case 
keeps a pole, the residue of which is slightly bigger than the perturbative one. 

As we shall see, the chiral solution will continue to possess these properties, even after 
we get a consistent equation for the quark propagator by the introduction of the transverse 
part of the vertex. 

11 



2. A multiplicatively renormalised equation 

The main effect of tlie inclusion of the transverse vertex is to shift the divergence from the 
integral of Ai to that of A2-F in Eq. (^4.3|) . Indeed, the integral over k^ of A2 in Eq. ( |4.1j ) 
now diverges logarithmically since A2{k^,q'^) behaves like 1/fc^ as k"^ — * 00, whereas the 
integral of the Ai term is finite. 



This fact means that the renormalised coupling 
by a finite renormalisation: 



is related to the usual one 



"/ 



asifif) 



l-as{^if)CF/A7TJdk^Zik^)A,{k^,^l}) 



(4.6) 



Hence we can now write our equation in terms of the true QCD coupling: 
1 



J^bM 



q'Fniq') 

Cf (ysifj'f) 



dk' 



+ Zn{k')A,{e,q')Fn{e + q'] 



Att 

„2\ A /1„2 „2\ T? /i„2 



(4.7) 



Setting g^ = /i^ in the above, we obtain: 



Cf asifJ'f) 



^bM 



1 - z^^iiiiLii^ I de 
Atc 



ZR{k^) 



Al(A;^/i2) 



^l) 



Znik')A2ie,fij)FRik' + fij) 



We can then equate both expressions for 1/J^h{fif) to obtain: 
1 = q'Fniq') 



Cf asifJ'f) 

A-K 

Cf asjl^f) 
47r 



dk^ 



2 ,,2\-\ 



Zn{k')A^{k\q')Fn{q')-Zn{k' 



Ai(P,/i^ 



2 ,.2^ 



dk'Znik') A2ie,q')FR{e + q')-A2ik',fij)FR{e + fij) (4.9) 



As can be seen the large k"^ behaviour of A2 cancels out and the integral is now ultraviolet 
convergent. 



5. The chiral solution 



5.1. Asymptotic behaviour 



The first test of the consistency of our results is a direct comparison of equation (|4.g| ) 
with perturbation theory, i.e. for g^ -^ 00. As we have explained, before subtractions, the 
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equation is ultraviolet divergent because of the terms proportional to A2. After subtraction, 
the leading logg^ conies from those terms in A2 that behave like l/k"^. The equation then 
becomes, for g^ ^ 00: 



l^ q^Fniq')- 



Cp ots f^^ 3dk 



All 



2P 



-FR{k'') Znik 



(5.1) 



Writing q^Fniq^) ^ 1 + ^\og{q^), Zr 
^ ^ as, provided that 



e 



1, one then gets a consistent solution to order 
3Cf as 



Svr 



(5.2) 



in agreement with one-loop results (in axial gauge). 



Furthermore, the equation also agrees with RG-improved perturbation theory. Assum- 
ing that for large q"^ one has Z^lq'^) ^ log(g^)~^, one then gets a consistent asymptotic 
behaviour q^F^^q^) ~ C + log(g^)~^, with C a constant, and ^ still given by Eq. (|5.2| ). 
Hence we see that Eq. (|4^ ) encompasses our perturbative knowledge of propagators. 



Finally, as we shall now explain, it suggests that the G = solution cannot be valid for 
arbitrary values of asifi/)- 



5.2. Critical value of as and chiral symmetry breaking 

By shifting the k"^ integration in the terms proportional to A2, one can recast the equation 
in the following form: 



with: 



0(g') 

/c(g^ k') 



):Fn{q') = <P{fi}) + J dk'lC{q\ k^)Tn{k' 



1 - ^^^ I de 

A-K 
Cp as 1 
A-n fc2 



2,Ai(P,g2^ 



ZR{k 



Zn{k' - q')A,{k', q') - Zn{k' - ^^j)A,{k\ /xj) 



(5.3) 



(5.4) 



In the case where (/>(g^) 7^ for all g^, the equation can be reduced to a Fredholm 
equation of the second kind. Its kernel /C is integrable and bounded, so that there is a 
unique solution. We obtain it by discretising Eq. (|5.3D and inverting the matrix equation 



thus obtained, which is equivalent to Fredholm's solution [0. We then get a smoother 
solution by introducing the obtained points as the input of an iterative method, where we 
use the left-hand side of Eq. ( p.3|) as the output, and the left-hand side as the input, which 

^ for all 



is the Liouville-Neumann method [|18l. This converges nicely as long as 
q'- 
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Figure 3: Massless quark solutions to the Dyson-Schwinger equations exist in the shaded region, 
the boundary of which is the value of as for which the equation becomes singular for some value 
of q^ , given the mass squared a which enters the Kallen-Lehmann representation of the gluon 
propagator. The thick curve shows the values of as at which the effective quark coupling ( [4 .21 ) 
becomes infinite. The curves are for r/ = 1. 



When 0(g^) has a zero, we have a Fredholm equation of the third kind. In that case, 
both the Fredholm solution, and the Liouville-Neumann iterations fail. This is because, 
effectively, the Fredholm solution involves integrals of J dq'^K{q'^,k'^)/(j){q'^), which are ill- 
defined, whereas the Liouville-Neumann solution involves a series with terms of the form 
J dk'^K{q'^, fc^)0(/c^)/(/)(g^), which clearly diverges near the zero of 0. 

In fact it is clear that the solution dramatically changes. When (pifi'j) = 0, the solu- 
tion J^ can at best determined up to a constant as the equation becomes homogeneous. 
Furthermore, it is not in general possible to remove the pole that the zero of (p introduces 
in JF, hence the propagator develops an imaginary part, which is not allowed in the t- 
channel. 0(g^) is 1 at a^ = and steadily decreases until it reaches a zero value (Note that 
Ai(fc^, g^) — i> as g^ — i> 0, hence is finite for all g^). becomes zero at small g^ first, and 
once 0(0) has crossed zero, then there will be a zero of 0(g^) at some nonzero g^ for larger 
values of as- Therefore there is a critical value of as past which is not positive definite. 

Physically, one can view the divergence of the kernel as the divergence of the effective 
coupling Q;/(g^) of Eq. ( [4.2| ) at some value of g^. This divergence suggests that the chiral 
equation stops having physically relevant solutions, and hence that chirality needs to be 
broken past a certain critical value of the coupling. One should really speak of the value 
of Q;5(/i)Z/j(yu), as one can always change the value of as according to Eq. ( |2.6| ). This 
interpretation is reinforced when one realizes that the singularity of the kernel is entirely due 
to the fact that we have neglected the terms proportional to G. If these are reintroduced, 
the singularity will disappear, and one will keep a solution. Hence, the singularity is the 
place at which the quarks develop a mass. 

We can in fact study this for a general gluon propagator, and hence not limit ourselves 
to the solution found in ref. |I[]. We simply need to assume that the gluon propagator has 
a Kallen-Lehmann representation: 

^^1^^ (5.5) 

g"^ J q"^ + a 

4>{q'^) can then be written, using the condition Z[iig) = 1: 



0(g') = 


j dap{a)^{a,q^) 


with $((T, g^) = 


\ f^l Cf as 



j deA,{k\q-)^\ (5.6) 



k^ + a ^ 

We show in Fig. 3 the region in the {as-, o'/fi'^) plane in which with $(o", g^) is positive 
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Figure 4: The solutions for massless quark propagators that come from the nonperturbative gluon 
propagator of ref. [0, for as{mf) = 0.2 (plain), 0.6 (dashed), 1.0 (dot-dashed) and 1.4 (dashed). 

for all g^. This means that if p(cr) has support in an interval [(To,o"i], then there will be 
massless solutions if as < a5"**™'(cri). 

We have plotted Fig. 3 for the function r] of Eq. ( p.l2|) equal to one. It is easy to see 
that a critical as will exist for a broad range of n and that its value can be calculated from 
Fig. 3. Indeed, as shown in the Appendix, for rj f 0, , 2^^2-)2 ) , we can write Ai = r^Af + Af. 
Unless t] is such that the sign of 



.,|Af(fcV<(=) + ^Af(fc7,')] 



^""-i^'' V^+V ,5.7) 

changes, there will exist a critical as- Its value will be adr]) = adrj = l)I{ri = 1)/I{ri). 
We therefore conclude that there exists a wide range of choices for the transverse vertex 
which give rise to chiral symmetry breaking. 

It is of course true that even for $ negative in part of [ctq, ai] it is possible to get > 0, 
hence the condition is sufficient only: in the shaded region, there will be massless solutions. 
It is interesting to note that large-cr modes in the Kallen-Lehmann density of the gluon 
imply chiral symmetry breaking in the quark sector. Only for a specific gluon propagator 
can one find the exact value of as beyond which these solutions do not exist anymore. This 
is what we are going to do in the next section. 



5.3. Chiral solution for a specific gluon propagator 

As explained in Section 2, the gluon propagator (|2.5| ) that we found in ref. 0] has the 
problems linked with the breakdown of multiplicative renormalisability. The procedure 
used here to obtain a consistent ansatz for the transverse part of the vertex can be extended 
to the gluon case. However, one may hope that the effect on gluons will be less dramatic 
than that on quarks. Indeed, we found that imposing the true asymptotic behaviour as 
g^ — * oo does not appreciably change the behaviour of the propagator at moderate q^ 
(compare the solution of the first paper of ref. 1 with Eq. ( p.5|) ). 

Hence we shall use that solution in the quark equation as an example of what the 
chiral solution looks like. We shall assume that both equations are renormalised at the 
same point, and that the two couplings have the same value. As we already explained, 
one should be able to derive a gluon propagator for any value of as{^J) as it can always 
be changed according to Eq. ( |2.5|) . As our gluonic equation did not respect multiplicative 
renormalisability, we were not able to do so, and got an optimum as ~ 1.4. To see what 
the effect of the coupling on the quark propagator is, we shall vary as independently of Z, 
although the two are really correlated. 

We show the result of this exercise in Fig. 4. As the value of the coupling grows, 
the quark propagator is enhanced near the origin. As as grows further, the propagator 
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/„2_|_ 2\2 

Figure 5: Same as Fig. 4, but for with a transverse vertex multiplied by ^44 . 

experiences oscillations, until the chiral solution is lost. It is interesting to note that the 
value of as that we obtained for the gluon propagator is very close to the critical value 
beyond which the chiral solution disappears. We also see that the criterion of Section 5.2 
can in practice work backwards, i.e. it not only predicts when there will be chiral solutions, 
but also when these will disappear. 

In order to investigate the sensitivity of the solution to our choice of the function D 
in Eq. ( 3.11|) , we consider two possible choices for that function. In Fig. 5, we show the 



behaviour of our solution for for V = q^ + p'^, instead oi V = (g^ + p'^Y (as in Fig. 3). As 
expected, the two solutions have the same high-g^ behaviour and are equal at g^ = /i^- We 
see that the low-momentum behaviours vary by a modest amount, and that the solutions 
only differ in the region g^ ^ fi'j, and only when the integrand oscillates a lot. This is 
where we do not trust our approximation for the angular integrals and hence the result is 
stable where our approximations hold. 



Figure 6: The dashed curve shows the value of Oc for various functions rj (see text), and the plain 
curve shows the value of the intercept of the propagator for as = 0.75 ac- 

Finally, Fig. 6 illustrates the insensitivity of our results to the choice of transverse vertex. 
We consider the functions rj = [{q'^ + p^)/{q^ + P^)^]" and plot the results in terms of n. We 
show the variation of the critical value of as-, which changes by a factor 2 when n changes 
by a factor 100. We also show the value of the intercept at the origin of the propagator 
calculated at as = O.TSctc, and which hardly changes with n. This clearly illustrates that 
the results we have obtained hold for a wide class of transverse vertices. 



6. Conclusion 



We have shown that the transverse part of the vertex plays an essential role in the quark 
DS equation in the axial gauge. It restores multiplicative renormalisability, allows chiral 
symmetry breaking and provides solutions which match with perturbation theory at large 
q^. Hence we have obtained an equation that possesses all the properties that are required 
in QCD. 

The solution that we obtain in the chiral phase has a pole at g^ = despite the fact 
that we expect quarks to be confined. This contrasts with the results of a similar analysis 
carried out in ref. ||15[ for the gluon propagator of ref. |^ in which it was shown that for the 
confining propagator of ref. [§] the quark propagator does not have a singular behaviour at 
g^ = and can therefore be considered to be confined. The present formalism can be used 
to solve for the chirality breaking solution, as we shall explain in a future publication [|16 



It is an interesting question to see whether these solutions are confined, and hence whether 
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confinement and chiral symmetry breaking are related. 

It has not so far been possible to obtain a solution to the (modified) DS equations for 
the gluon propagator for time-like momenta. This is because we expect that in this regime 
the propagator develops an imaginary part with cuts corresponding to the thresholds for 
glueball production etc. and consequently such a solution involves coupled non-linear equa- 
tions for the real and imaginary parts of ZRik!^). Therefore the corresponding solution for 
the quark propagator we have obtained is only valid for quarks propagating with space-like 
momenta (in the t-channel). It is tempting to try to analytically continue the solution ob- 
tained into time-like momenta by fitting the solution obtained to known analytic functions 
and possibly then using the DS equation for the quark propagator to extract information 
about the gluon propagator with time-like momenta. Unfortunately this has not been pos- 
sible. We have obtained several very accurate fits (errors nowhere worse than 2% ) to the 
curve shown in Fig. 4, using different parametrisations inside various analytic functions. 
Whereas these all fit the quark propagator remarkably well in the fit region they give wildly 
different projected behaviours for the quark propagator for time-like momenta. We there- 
fore have to accept for the moment that the important problem of describing the exchange 
of soft gluons and quarks in the s-channel remains unsolved. 

As we explained earlier, the equation for the gluon propagator also will have to be 
modified to incorporate a transverse vertex and recover multiplicative renormalisability. 



Furthermore, as has been pointed out in ref. |19] the DS equations for the quark and gluon 



propagators are coupled and after having obtained the solution for the quark propagator 
one must check the assumption that one may neglect quark loops in the equation for the 



gluon propagator. In ref. []19| it was shown that for the gluon propagator of ref. [§ this was 
not the case and the gluon propagator had to be modified accordingly. A similar analysis 
should also be carried out for the gluon and quark propagators discussed here. 
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7. Appendix: Kernels 



The kernels entering Eq. ( ^.l] ) are given by the following expressions, where we define 



For the BZ equation (p.6|): 

Ai(p) 



|_,„_4)^,_i(^4 



and 



A2(P) 



;i-p)^+%-i)(^(i + V^' 



+^(p-4)(l + p)^l-i(f + l 



(7.1) 



(7.2) 



n2„2 



For Eq. ( p.l5| ), and for r/ a function of ^f^l^y ~ (2^^ only: 

p2 + 16 + 2p' 



Ai(p) 



ri_ 
Ap 
1 

+ 4 



-P-8 + 

-p' + 2 p + 
Vp'-4p 



v/;^^T4 



+ P 



^ ^<-^)i^^^^ 



(7.3) 



A2(p) = ^|9p + p2 + 8--^=^=[p3 + p2 + 18p + 16] 



+ 0{p-l) 

+ ^(1-p) 
+ e{A-p) 



2p 



-3 + p' + p 
2 

y + p^ + iop + 8 

4v'p^-4p 
Q — p'^ — p p^ — p^ — lOp - 



4Vp'-4p 



(7.4) 
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Abstract 

Within the framework of the Dyson-Schwinger equations in the axial gauge, we study 
the effect that non-perturbative glue has on the quark propagator. We show that Ward- 
Takahashi identities, combined with the requirement of matching perturbative QCD at 
high momentum transfer, guarantee the multiplicative renormalisability of the answer. 
Technically, the matching with perturbation theory is accomplished by the introduction of 
a transverse part to the quark-gluon vertex. We show that this transverse vertex is crucial 
for chiral symmetry breaking, and that massless solutions exist below a critical value of 
as- Using the gluon propagator that we previously obtained, we obtain small corrections 
to the quark propagator, which keeps a pole at the origin in the chiral phase. 
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1. Introduction 



The Dyson-Schwinger (DS) equations of motion are one of the main tools for the investiga- 
tion of non-perturbative effects. These equations are particularly suited to the study of the 
interface between perturbative and nonperturbative regimes, as they allow us to compute 
the evolution of Green functions when one enters the infrared region, and hence suggest 
modifications to the usual Feynman rules. 

In a previous paper, we have studied [I] the DS equations in the quenched approxima- 
tion, i.e. for pure gauge QCD, and found that, as was first suggested by Cornwall [2], some 
solutions for the gluon propagator are fiatter than a pole in the infrared region. In order 
to avoid the question of ghost propagators, we worked in the axial gauge. Other groups 
[3, 4, 5] have also found that fiatter solutions exist, in covariant gauges, although it is at 
present difficult to relate our results to theirs. The absence of a pole in the gluon propaga- 
tor at A;^ = is natural if one assumes that gluons do not propagate to infinity, i.e. these 
solutions should describe confined gluons. Furthermore, as was pointed out by Landshoff 
and Nachtmann [6], the existence of such solutions is highly desirable in phenomenological 
applications, and gives us a practical way to extend usual perturbative estimates to the 
strongly-interacting sector of the theory. The use of these solutions has already met with 
some success in diffractive calculations [7]. 

The DS equations constitute an infinite tower of integral relations between n-point 
functions. Hence, by themselves, they cannot be solved. For instance, if we consider the 
DS equation for the quark two-point function (the propagator), it involves a three-point 
function (the quark-gluon vertex) as well as the gluon two-point function (the propagator). 
One can imagine solving the pure-gauge DS equations for the gluon propagator as a first 
approximation, hence we shall assume here that the gluon propagator is known. However, 
resorting to higher-order DS equations to calculate the vertex will only bring in higher- 
order n-point functions, which are also a priori unknown. Hence one needs to use another 
input. 

The simplest assumption is to take the perturbative 7^ vertex, which is the so-called 
"rainbow" or "ladder" approximation. One can be however more sophisticated, and follow 
the observation of Baker, Ball and Zachariasen (BBZ) [8], namely that one needs to choose 
a vertex that will obey the Ward- Takahashi-Slavnov- Taylor identities. These constitute 
another nonperturbative statement of field theory, and hence must be valid in general. 

In our study of the gluon propagator [I], we chose for the three-gluon vertex the simplest 
function that would obey the Ward-Slavnov- Taylor identities. We showed that besides the 
original l/k"^ solution, there exists another solution, behaving like Ijk^'^ near k^ = 0. To 
obtain a full picture of nonperturbative effects in QCD, one then needs to calculate the 
quark propagator corresponding to that gluon propagator. As we shall see, the problem of 
multiple solutions arises again, as for sufficiently small 05, a chiral solution exists together 
with a massive one. We shall show that one can constrain the spectral density of the gluon 
propagator and the value of the coupling so that only the massive solution survives. 



As the equation for the quarks is much simpler than that for the gluons, we investigate 
in more detail the effect of the ansatz for the quark-gluon vertex on the quark propagator. 
Namely, the minimal vertex, used by BBZ for gluons and by Ball and Zachariasen (BZ) for 
quarks [f2], leads to the breakdown of multiplicative renormalisability. We demonstrate 
that, following the method of Curtis and Pennington (CP) [9], it is possible to recover 
it by choosing a specific form for the quark-gluon vertex, which agrees with perturbation 
theory in the ultraviolet region, and which generalizes the B(B)Z ansatz. We also show 
that massive solutions exist only for the CP vertex. Simpler ansatze lead to an inconsistent 
ultraviolet behaviour, unless the propagator is massless. We also explain under which 
conditions massless solutions will exist, even for the CP vertex. We proceed to solve the 
DS equation in the massless phase, and show that nonperturbative effects do not remove 
the pole of the quark propagator in the chiral phase of the QCD vacuum, suggesting that 
confined quarks can exist only after chiral symmetry breaking. 

This paper is set out as follows. In Section 2, we review our results for the gluon 
propagator, and outline the formalism that lead to these. We point out the problems linked 
with the BBZ choice of a three-gluon vertex. In Section 3, we discuss the DS equation for 
the quark propagator in the axial gauge, and the approximations made when imposing the 
Ward-Takahashi identity. In Section 4, we consider the part of the propagator that preserves 
chiral symmetry. We show that the problem can be reduced to a one-dimensional integral 
equation, which we then renormalise. We show explicitly that we recover multiplicative 
renormalisability. In Section 5, we give a criterion for the existence of massless solutions, 
and then proceed to the numerical solution of the equation in the massless case, for the 
gluon propagator that we previously derived. 



2. The Gluon Propagator 

The fundamental observation [8] is that one can use the Ward-Slavnov- Taylor (WST) iden- 
tities to obtain a closed equation from the DS equation, if one neglects quark loops and 
works with gluons only. In the axial gauge [n.A = 0), where we do not have to worry 
about ghost degrees of freedom, the DS equation relates the propagator to the three- and 
four-point vertices. 

The three-point vertex, FJf^ can be split into a part P^ transverse to the external 

gluon momentum k^, Tj ■ k = 0, and a longitudinal part, P^^ , with Tj^ ■ k ^ 0. As 
the vertex is a third rank tensor, depending on the three vectors entering the vertex, this 
split is not uniquely defined. The WST identities relate P^^ to the propagator, and these 
identities can be "solved", hence producing an ansatz for the longitudinal part [8, 12]. It 
is this ansatz that we define as being the longitudinal part of the vertex in the following. 

One can then obtain a closed equation for the propagator if one makes two assumptions. 
First, the full propagator is supposed to have the same spin and gauge structure as the free 
propagator: 



D,.{k') = -DliKn)^ (2.1) 

where 

DliK n) = V - ^'""^ ^!^'^' + ^^/^ (2.2) 

^ n.k [n.k)'^ 

This spin structure projects out the four-point vertex contribution to the equation. If one 
further assumes that the three-point vertex is dominated by its longitudinal part F^ , which 
is itself known from the propagator via the WST identities, one obtains an equation which 
involves only the gluon propagator. 

This equation still needs to be renormalised. One can show without approximation [8] 
that a propagator with the spin structure (2.1) has to be singular as k^ — )• 0. Thus the 
inverse propagator should vanish at the origin. This allows one to subtract the quadratic 
divergences, thus renormalising the gluon mass to zero and cancelling the tadpole graphs. 
One is then left with logarithmic divergences, which can be dealt with through wave func- 
tion renormalisation at a renormalisation point /i^, by defining 

Z{e) = Z{ij'^)Zn{e) (2.3) 

However, the resulting definition of the renormalised coupling constant ag is slightly dif- 
ferent from the usual one: 

where «& is the bare coupling and ^ a function linear in Z, the explicit form of which 
can be found in ref. [11]. It has to be noted at this point that Eq. (2.4) leads to two 
related problems. First of all, the propagator asymptotically behaves like 1/ log(A;^)^^'^^ 
as A;^ — )• oo [8], which is not what is obtained in perturbation theory. Furthermore, one 
obviously loses multiplicative renormalisability. This has the nasty consequence that, if 
the same truncation is applied to the quark propagator, the renormalised couplings are not 
equal anymore, because the functions ^ are different for quarks and gluons. 

The solution to the resulting equation, found in ref. [8], has a very singular behaviour 
proportional to Ijk'^ as the momentum squared goes to zero. This was advertised as a signal 
for a confining potential. It has to be pointed out that in order to renormalise the equation 
in the presence of a Ijk'^ singularity, one needs to perform additional subtractions, and 
hence one obtains such solutions only if one assumes they exist. In ref. [1], we exploited 
the fact that the integral equation obtained in ref. [8] is non-linear and can in general 
admit several solutions. We solved the original equation, without additional subtractions, 
and obtained a gluon propagator which, although singular as A;^ — )• as required in the 
axial gauge, has a cut singularity, fiatter than a pole. This describes a soft gluon which 
is confined rather than confining and can conveniently be used in models of the pomeron 
such as that suggested by Landshoff and Nachtmann [6]. This propagator may be written 



in the axial gauge [n.A = 0) as 



Znik' 



h.2 



\^A^" 



+ c In 






+ e 



(2.5) 



The dimensionless constants ai, 61, 02, ^2, c, d and e, are given in the Table. The mass scale 
fig is the value of momentum at which Zji^k^) takes the value 1, so that the propagator 
matches the value of the free gluon propagator. It is not obtained from the DS equation, 
which is scale invariant, but rather from a fit to the total cross-section for proton-proton 
scattering calculated using this propagator. An optimum value of 0.8 GeV was found for 
/ip. It is worth pointing out that most of the ideas discussed here can be applied to other 
theories, living at other mass scales, such as technicolour, and that one would only need to 
adjust /ip to do so. 

Table: The constants used in Eq. (2.5) 



ai 


0.88 


c 


0.59 


02 


-0.95 


d 


2.1 


w 


0.22 


e 


4.1 


h 


0.86 


agiflg) 


1.4 



Note that in the case of a multiplicatively renormalisable theory, the actual value of as 
should not matter. In principle, one could always rescale as by varying fi according to 



as{fi')ZR{fi''^,q'^) 



as{p)ZR{p^,q^) 
Zn{fi\fi")ZR{fi'\q') 



(2.6) 



where the first argument of Z refers to the renormalisation scale: Z[fi'^^ fi^) = 1. However, 
our gluonic solution does not obey these equations because of (2.4), and one obtains a 
solution only for a given value^. This value of ag[fig) enters phenomenological calculations, 
which then determine fig. 

Hence the assumption of a given spin structure (2.1) and the neglect of the transverse 
part of the three gluon vertex leads to a DS equation which admits two solutions, and 
which breaks multiplicative renormalisability. In order to explain the discrepancy of the 
BBZ result with more recent ones, the assumed spin structure of the gluon propagator has 
been put into question [5]. However, there is no doubt that such solutions exist. Further 
solutions can of course exist, given that the equation is highly nonlinear. The neglect of 
the transverse part of the vertex, on the other hand, may have important consequences, 
which we shall explore in this paper for the quark propagator. 



^Note that there was a mistake of a factor 2 in our definition of ag in Refs. [1] 



3. The Dyson- Schwinger equation for quarks 

Following their work with Baker on the gluon propagator [8], Ball and Zachariasen have 
considered the DS equation for the quark propagator [12]. Choosing again a purely longi- 
tudinal quark-gluon vertex, they found that the quark propagator corresponding to a Ijk'^ 
gluon was suppressed, and became constant near the origin. Furthermore, they showed 
that a chiral solution was always possible, and argued that a solution could be found 
which would break chiral symmetry. We shall see that, in the case of a less singular gluon 
propagator, only the chiral solution is possible if one neglects the transverse part of the 
vertex. However, implementing the improvements proposed by Curtis and Pennington [9], 
we shall show that multiplicative renormalisability (and hence the transverse part of the 
vertex) leads to an equation that allows chiral symmetry breaking, and even suggests a 
range of parameters (describing the coupling and the gluon propagator) which would lead 
to a unique, massive, solution. 



3.1. The Ball- Zachariasen equation 

We define the propagator for a quark with momentum q as: 

S{q) = F{q')^-q + G{q') (3.1) 

so that F[q^) represents the chiral-symmetry conserving part and G{q^) represents the 
chiral-symmetry breaking part. We shall also use the equivalent notation: 

SW) = ,2_s2(,2) (3-2) 

The free propagator is therefore obtained by setting F[q^) to l/q^ and G{q^) to zero 
(neglecting all current masses). In the axial gauge Fq. (3.1) is not the most general form for 
the quark propagator. First of all, the functions F and G can depend on n.q as well as q^. 
As in the case of the gluon propagator [8, 1], we seek solutions for which these functions are 
independent of n.q. Furthermore, the gauge dependence can also arise through extra spin 
structures, proportional to 7 • n and n.q^ which although absent for the free propagator can 
in general occur for the dressed propagator. As pointed out in ref. [12]), these extra terms 
drop out of the equations for F and G if one specializes to a gauge vector orthogonal to q. 
This gauge choice n.q = is the only one that makes the algebra sufficiently tractable that 
reliable solutions to the DS equations can be obtained. The question of the dependence of 
the quark propagator on n, although important, is difficult to address in the axial gauge, 
and beyond the scope of this paper. 

The DS equation for S{q)^ sketched in Fig. 1, in Fuclidean space is given by 

1 = 7 • (? Sb{q) - Cpa^j -^^Z,{e)rD,,{k, n)Sb{q - k)T''{q - k, q)Sb{q) (3.3) 



k=q-p 




Figure 1: A pictorial representation of Eq. (3.3). The hatched circles represent the exact two-point 
functions, and the cross-hatched circle the exact three-point function. 

where the subscript b on a^ Z and S indicate that these are bare quantities which will have 
to be renormalised. Cf = 4/3 is the quark Casimir invariant. 

In general it is not possible to solve this equation because it involves the unknown quark- 
gluon vertex function T'^[q — k^q). This vertex function is related to a four-point Green 
function via a DS equation. Thus we get an infinite tower of coupled integral equations. 
On the other hand, in the axial gauge, this quark-gluon vertex function connecting a quark 
of momentum p = q — k to one of momentum q through the absorption of a gluon of 
momentum A;, obeys the Ward-Takahashi (WT) identity sketched in Fig. 2: 

k^T^{p,q) = S{q)-' - S{p)-' (3.4) 




p 



Figure 2: A representation of the Ward-Takahashi identity of Eq. (3.4) using the same convention 
as for Fig. 1. 

One can then solve this equation to determine F^, the part of the vertex longitudinal 
to k. Following ref. [12], we obtain: 



S{q)T^(p,q)S{p) = -[F,{p) + F,{q)]Y + 7:mp)-F,{q)] 



2 ' '"' '^'' ' 2' '"' '^'' p2 _q2 

+ lmp)-F,{q)]^4^^' 

2 p'^ — q'^ 

+ [G,{p)-CMq)] ^'^'l^\'^^' (3.5) 

pz — q'^ 

The simplest approach consists, as suggested in ref. [13], in neglecting the transverse part 
of the vertex: one then assumes that the integral equation is dominated by the part (3.5) 
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which is determined from the above WT identity.^ The resulting integral equations for 
the two functions Fb[q'^) and Gh{q^) separate and we have 



1 = q Fb{q )-Cf ab I j^^^Zb{k )D^^{k,n)< ^ 



jJ.U 



+ 
and 



{q' 



^ [Fb{q') - Fb{p')] (^hn,, + q,p, + p.q^i \ (3.6) 



Gb{q') -Crab I -^^Zb{P)D,,{k,n) 



[q^ — p^) 



Gb{q ) - Gb{p )\{q- kS^i, + q^pi, + p^q^,) ) (3.7) 



where again we have introduced the subscript b on F and G to indicate that these quantities 
are to be renormalised. (There is no renormalisation for E since this function is zero for 
the free propagator, i.e. in the absence of an explicit current mass there is no parameter 
in the QCD Lagrangian to renormalise, hence the ratio F/G does not get renormalised). 



3.2. Multiplicative renormalisability and the transverse part of 
the vertex 



As was observed in [12], Eq. (3.6, 3.7) are not multiplicatively renormalisable. In other 

4) 



words, if we renormalise the quark and gluon wavefunctions by imposing Fb{q^) = Fb{fi'j 



Fniq^), Gb{q^) = Fb{jJ^f) GR{q^) and Zb{q^) = Zb{jJ^g)ZR{q^), additional terms need to be 
introduced in the definition of the renormalised coupling, as in Eq. (2.4). These terms 
are not the same in the quark and in the gluon case, so one loses not only multiplicative 
renormalisability, but also the universality of the QCD coupling constant! 

These problems come from the ultraviolet region, and can be traced back to the neglect 
of the transverse part of the vertex in the solution of the WT identities: the vertex (3.5) 
does not match the perturbative one at high momentum transfer. Hence, it is necessary 
to postulate a transverse part that will restore multiplicative renormalisability. Curtis 
and Pennington have shown [9] that this goal can be achieved in QED by considering the 
perturbative limit of the vertex. Their argument can be trivially extended to QCD in the 
axial gauge: the one-loop corrections to the propagator are identical to those of QED up to 
the quark Casimir invariant Cp, and the vertex corrections are also the same: because the 
vertex and the wavefunction renormalisation constants Zi = Z^ are equal in this gauge, 
and because Z^ is a function of Cp only, the diagrams involving the three-gluon vertex. 



^In other gauges the Ward identity Eq. (3.4) must be replaced by the full Slavnov- Taylor identities [14] 
which involve the interactions of Faddeev-Popov ghosts. Nevertheless it was argued in ref. [15] that the 
contributions from Faddeev-Popov ghosts are small and so the Ward identity of Eq. (3.4) was imposed on 
the integral equation formulated in a covariant gauge. 
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which depend on Cf — CaI'^, with Ca the adjoint Casimir invariant, have to canceL Hence 
the argument is totally similar to that of ref. [9] and goes as follows. 

In general, the quark propagator has the perturbative limit 

lim :F{q^) = l + '^\n^ (3.8) 

with ^ the anomalous dimension. The one-loop vertex can be shown [9] to tend to 



lim Vl-\p,q) = 7, 



1 «^^i 1^"" 



- 1 i — \1 -Villi -V- (liil -V^ifliil -V) (3.9) 

47r q-^ 

One can compare this answer with the vertex (3.5) in the same limit, to conclude that 
the transverse part must give: 



2 lim r (p,5) = — - — - — -{-qi.i-v^i-Tii.i-v) (3.10) 

This tensor is indeed transverse to q ^ k. One can then extend the tensor structure of 
this vertex so that it becomes transverse to k for any value of p and q. We find that the 
simplest extension is: 

c^( ^rT( wr^ (T( ^ ^i ,, 7.(g^-p^)-(g.+P.)(7-g-7-p) ,,.., 

SW^ [P, q)S{q) = {J^{q) - J^ip)) (3.11) 

where X) is an expression symmetric in p and q, and behaving as q^ in the Xaxge-q^ limit. 
X) must not introduce any singularity, must be symmetric in q and p, an must satisfy 
T^{p,p) = 0. In general, 

\{q^+p^) [q^+p^Y) 

with J] a regular function such that ?y(0,0) = 1. 

It is worth noting that the tensor structure of the transverse vertex (3.11) can be 
generalised. Indeed, Ball and Chiu [10] have given a set of 8 independent tensor structures 
spanning the space of regular transverse vertices. Only three of these have the correct 
helicity structure to contribute to the DS equation in the chiral limit, and are (in the 
notation of Ball and Chiu): 

T^ = [p^{k.q)-q^{p.k)h-{p + q) 

T^ = k^Y-k^l-k 

n = r{/-p')-{q + Pri-{q-p) (3.13) 
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The transverse vertex given by Eq. (3.11) corresponds to: 

where the eUipses refer to terms that vanish in the DS equation. This differs from the 
structure of ref. [9], where the transverse vertex was chosen to be proportional to Tq. This 
can be understood from the fact that the expression of the longitudinal part of the vertex 
used here differs from theirs, by a transverse tensor proportional to Ts, in such a way that 
we get the same high-(j'^ or p^ leading term for the total vertex. The addition of such extra 
subleading terms is always possible, and we shall investigate the effect of these through the 
variation of function T> in Eq. (3.11). 

The transverse vertex (3.11), together with the longitudinal one (3.5), can then be 
included in Eq. (3.3) to obtain: 



1 = p'F,ip) + CFa,Jj-^DlZ,ie 



{p^Fb{p)- q^Fb{q)){6f,^{q^ -p^)+Pf,p^- qf,q^+Pf,q^-p^qf,)\ 
+ ^ 1 (3.15) 

4. Solving the equation 
4.1. The angular integral 

The angular integral in Eqs. (3.6, 3.7, 3.15) cannot be performed analytically without 
prior knowledge of the function F, since this occurs with the argument [q — kY . The 
approximation, first proposed by Schoenmaker [11] in the gluon case, consists in replacing 
a function f{{q — kY) by f{q^ + k"^). This is clearly valid in the regions k^ » q^ and 
k^ « q^. One then needs to find a function / which is sufficiently slowly varying that the 
error generated by the integration over the whole range of k^ be negligible. 

The first choice would he f = F . In this case, the contribution of the longitudinal part 
of the quark-gluon vertex to the DS equation vanishes. Hence the BZ equations (3.6, 3.7) 
admit only the trivial solution F[q^) = l/q^^ G = 0. If the transverse part of the vertex 
is added, the resulting equation (3.15) leads to a propagator close to the perturbative one, 
with a pole at k^ = 0. Neither solution is a slowly varying function, hence the approximation 
is not justified in this case. 

In the following, we shall use j = T . Our solutions show that this assumption is 
reasonable so that the above ansatz is justified a posteriori. The only region where this 
approximation may have led to substantial errors would be for small q^ sa A;^, but then 
the integrand is suppressed because (as can be seen from Eq. (2.5)) the function Zji^k^) 
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vanishes. This approximation keeps the contributions both from the longitudinal and 
from the transverse part of the vertex, hence we shall obtain nontrivial solutions both 
for Eqs. (3.6, 3.7) and (3.15). 

With this approximation the angular part of the integration over momentum k in 
Eqs. (3.6, 3.7, 3.15) may be performed, and both equations can be recast into the fol- 
lowing form if G = 0: 

1 = q'F,{q') 

47r 



dk' Z{e)A^{e/q') F,{q') + Z{k')A,{k'lq')F,{k' + q')\ (4.1) 



The kernels Ai and A2 are given, both in the case of Eqs. (3.6, 3.7) and in the case of 
Eq. (3.15) in the Appendix. This is a Eredholm equation of the second kind which must 
now be solved numerically. Before we can do this, however, we must consider the question 
of the renormalisations required to absorb the ultraviolet divergences. 



4.2. Renormalisation 

The renormalisation procedure is similar whether one neglects the transverse part of the 
vertex or not. The renormalisation constant Zb[fi'^) has been introduced in Eq. (2.3). 
As mentioned in Section 2 we choose this constant such that the renormalised function, 
Zji^q^) takes the value 1 ai q^ = ij\ with //^ taken to be 0.8 GeV. We likewise introduce a 
renormalisation constant J-'b[fif) such that 



and 

Gr{/ 



G,{q') 



are ultraviolet finite and we choose it so that the quantity J^R^q^) takes the value 1 at 
q^ = fi^f. Eq. (4.1) then becomes finite if one rewrites it in terms of the coupling 



1 - abZ{iif)CF/^Tr I dk^Z{k^)Ai{k\ l^j) 



J 1 „ o'/,, \ry /A- r J1.-2 crti.-2\ \ ii„2 ,,2\ \^-^) 



and the renormalised equation (4.1) becomes: 

1 = {l-^^jdk'Zn{k')[A,{k\q')-A,{k\,j',)\)Fn{q') 

^^ "^- j dk'Znik') (A2(^^ q')Fnik^ + q') - A,{k\ ^})Fnik^ + ^})) (4.3) 



47r 
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1. Renormalised BZ equation 

In the case of the BZ kernels, the UV divergence of the equation is concentrated in the Ai 
term, which has a log singularity as k^ — )• oo. This means that the coupling (4.2) cannot 
be related to the usual renormalised QCD coupling 

asifif) = Zb{fif)ab (4.4) 

Hence one loses the universality of the QCD coupling. 

Furthermore, the equation for the chirality breaking term G is identical to Eq. (4.3) 
with the replacement F — )• G, 1 — > 0. The leading behaviour at high q^ comes from the 
remnant of the UV term. If we assume, in agreement with the renormalisation group, that 
^-r(<?^) ~ l/[9^ log(<?^)] ^^ <?^ ~^ °°7 the leading terms of the equation for the G term are: 

^ Gn{q') - %^ log{log{q'))Gn{q') (4.5) 

47r 

which has only G = as a consistent solution. Hence, unless the gluon propagator is highly 
singular, the BZ equation does not lead to chiral symmetry breaking. 

If one overlooks these problems, and goes ahead to solve the equation, one obtains a 
function Fji^q^) which is very close to the perturbative Ij q^ ■ The propagator in this case 
keeps a pole, the residue of which is slightly bigger than the perturbative one. 

As we shall see, the chiral solution will continue to possess these properties, even after 
we get a consistent equation for the quark propagator by the introduction of the transverse 
part of the vertex. 

2. A multiplicatively renormalised equation 

The main effect of the inclusion of the transverse vertex is to shift the divergence from the 
integral of Ai to that of A2-F in Eq. (4.3). Indeed, the integral over k^ of A2 in Eq. (4.1) 
now diverges logarithmically since ^2{k'^^(f') behaves like Ijk^ as k^ — )• 00, whereas the 
integral of the Ai term is finite. 

This fact means that the renormalised coupling (4.2) is related to the usual one (4.4) 
by a finite renormalisation: 

^ l-as{iif)CF/4:TrJdk^Z{k^)Ai{k\li}) ^'^ 

Hence we can now write our equation in terms of the true QCD coupling: 



^ - q'Fni/) 



Cf ctsil^f] 



Jdk'[Znik')A,ik\q')Fniq' 



+ Zn{e)A2{k\q')Fn{e + q')] (4.7) 
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Setting q^ = fij: in the above, we obtain: 



1^ = 1 - ^^LMtll I dk' 

J^bifif) 47r 






.8) 



We can then equate both expressions for l/J-'b[fif) to obtain: 

1 = q'Fniq') 

dk 



CpasM /,,/ 



47r 
47r 



Znie)A,ik\q') Fniq') - Znik')^'^^''^''^ 



fij 



I deZnik') [A2{k\ q')Fn{e + q') - A,{k\ fi'j)Fn{e + /i^)] (4.9) 



As can be seen the large k^ behaviour of A2 cancels out and the integral is now ultraviolet 
convergent . 



5. The chiral solution 



5.1. Asymptotic behaviour 

The first test of the consistency of our results is a direct comparison of equation (4.9) 
with perturbation theory, i.e. for q^ — )• 00. As we have explained, before subtractions, the 
equation is ultraviolet divergent because of the terms proportional to A2. After subtraction, 
the leading logq^ comes from those terms in A2 that behave like 1/k^. The equation then 
becomes, for (j'^ — )• 00: 

1 - I'Fniq') - %^ r^ ^-^Fnik-) IZniP)] (5.1) 



47r J^ 2A;2 

Writing q^Fji^q^) sa 1 + ^\og{q^), Zr ~ 1, one then gets a consistent solution to order 
^ ~ as, provided that 

( = ^ (5.2) 

in agreement with one-loop results (in axial gauge). 

Furthermore, the equation also agrees with RG-improved perturbation theory. Assum- 
ing that for large q^ one has Zji^q^) sa \og{q^)~^ , one then gets a consistent asymptotic 
behaviour q^FR{q^) ~ C -|- \og{q^)~^ ., with C a constant, and ^ still given by Eq. (5.2). 
Hence we see that Eq. (4.9) encompasses our perturbative knowledge of propagators. 

Einally, as we shall now explain, it suggests that the (7 = solution cannot be valid for 
arbitrary values of as{fif). 
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5.2. Critical value of as and chiral symmetry breaking 

By shifting the k^ integration in the terms proportional to A2, one can recast the equation 
in the following form: 



with: 



q'):FR{q') = <P{fi-j) + I dk^]C{q\ P)J^r{P) (5.3) 

Ai(P,?2) 



(/) = 1- ^jde 



Mk' 



q^ 



K{q\k') = ^^^[Znie-q')A,ik\q')-Znie-^})A,ik\^})] 



(5.4) 



In the case where (f>{q'^) 7^ for all q^^ the equation can be reduced to a Fredholm 
equation of the second kind. Its kernel /C is integrable and bounded, so that there is a 
unique solution. We obtain it by discretising Eq. (5.3) and inverting the matrix equation 
thus obtained, which is equivalent to Fredholm's solution [17]. We then get a smoother 
solution by introducing the obtained points as the input of an iterative method, where we 
use the left-hand side of Eq. (5.3) as the output, and the left-hand side as the input, which 
is the Liouville-Neumann method [18]. This converges nicely as long as (f>{q'^) 7^ for all 

When (f>{q'^) has a zero, we have a Eredholm equation of the third kind. In that case, 
both the Eredholm solution, and the Liouville-Neumann iterations fail. This is because, 
effectively, the Eredholm solution involves integrals of J dq^K[q^^k^)/(f>[q^)^ which are ill- 
defined, whereas the Liouville-Neumann solution involves a series with terms of the form 
J dk^K[q^^ k'^)cf){k'^) I cf){q^) , which clearly diverges near the zero of cf). 

In fact it is clear that the solution dramatically changes. When (f>{fi'j) = 0, the solu- 
tion J-' can at best determined up to a constant as the equation becomes homogeneous. 
Eurthermore, it is not in general possible to remove the pole that the zero of (f> introduces 
in JF, hence the propagator develops an imaginary part, which is not allowed in the t- 
channel. (f>{q'^) is 1 at 05 = and steadily decreases until it reaches a zero value (Note that 
Ai(A;^, q^) — )• as (j'^ — )• 0, hence (f> is finite for all q^). cf) becomes zero at small q^ first, and 
once ^(0) has crossed zero, then there will be a zero of cf>{q^) at some nonzero q^ for larger 
values of as- Therefore there is a critical value of as past which cf) is not positive definite. 

Physically, one can view the divergence of the kernel as the divergence of the effective 
coupling af{q^) of Eq. (4.2) at some value of q^ . This divergence suggests that the chiral 
equation stops having physically relevant solutions, and hence that chirality needs to be 
broken past a certain critical value of the coupling. One should really speak of the value 
of as{fi)ZR[fi)^ as one can always change the value of as according to Eq. (2.6). This 
interpretation is reinforced when one realizes that the singularity of the kernel is entirely due 
to the fact that we have neglected the terms proportional to G. If these are reintroduced, 
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10 10 10 10 

Figure 3: Massless quark solutions to the Dyson-Schwinger equations exist in the shaded region, 
the boundary of which is the value of as for which the equation becomes singular for some value 
of q^, given the mass squared a which enters the Kallen-Lehmann representation of the gluon 
propagator. The thick curve shows the values of as at which the effective quark coupling (4.2) 
becomes infinite. The curves are for r] = 1. 



the singularity will disappear, and one will keep a solution. Hence, the singularity is the 
place at which the quarks develop a mass. 

We can in fact study this for a general gluon propagator, and hence not limit ourselves 
to the solution found in ref. [1]. We simply need to assume that the gluon propagator has 
a Kallen-Lehmann representation: 



da 



q'^ J q^ -\- (T 

[q^) can then be written, using the condition Z[fig) = 1: 



(5.5) 



^{/) = 


j dap{a)^a,q^) 


with$(cr, 5^) = 


\ 1^1 Cpas 
Xlil + a 47r 



j^^^^i^.e,^ 



(5.6) 



We show in Fig. 3 the region in the (05, cr///^) plane in which with ^{a^q^) is positive 
for all q^ . This means that if /o(cr) has support in an interval [(To,cri], then there will be 
massless solutions if as < ^^'^'^"'(cri). 

We have plotted Fig. 3 for the function rj of Eq. (3.12) equal to one. It is easy to see 
that a critical as will exist for a broad range of rj and that its value can be calculated from 
Fig. 3. Indeed, as shown in the Appendix, for r] (O, fJJ'^2)2 ) , we can write Ai = r]A^ -\-A^. 
Unless T] is such that the sign of 



'' (g2+p2)2 



Hv) 



dk 



,A'^{P/q') + ^A^{P/q')] 



(5.7) 



(F + a) 

changes, there will exist a critical as- Its value will be ac{r]) = ac{r] = l)/(?7 = l)/I{r]). 
We therefore conclude that there exists a wide range of choices for the transverse vertex 
which give rise to chiral symmetry breaking. 
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Figure 4: The solutions for massless quark propagators that come from the nonperturbative gluon 
propagator of ref . [1], for as(mf) = 0.2 (plain), 0.6 (dashed), 1.0 (dot-dashed) and 1.4 (dashed). 

It is of course true that even for $ negative in part of [ctq, cti] it is possible to get ^ > 0, 
hence the condition is sufficient only: in the shaded region, there will be massless solutions. 
It is interesting to note that large-cr modes in the Kallen-Lehmann density of the gluon 
imply chiral symmetry breaking in the quark sector. Only for a specific gluon propagator 
can one find the exact value of as beyond which these solutions do not exist anymore. This 
is what we are going to do in the next section. 



5.3. Chiral solution for a specific gluon propagator 

As explained in Section 2, the gluon propagator (2.5) that we found in ref. [1] has the 
problems linked with the breakdown of multiplicative renormalisability. The procedure 
used here to obtain a consistent ansatz for the transverse part of the vertex can be extended 
to the gluon case. However, one may hope that the effect on gluons will be less dramatic 
than that on quarks. Indeed, we found that imposing the true asymptotic behaviour as 
(j'^ — )• oo does not appreciably change the behaviour of the propagator at moderate q^ 
(compare the solution of the first paper of ref. 1 with Eq. (2.5)). 

Hence we shall use that solution in the quark equation as an example of what the 
chiral solution looks like. We shall assume that both equations are renormalised at the 
same point, and that the two couplings have the same value. As we already explained, 
one should be able to derive a gluon propagator for any value of as{fi) as it can always 
be changed according to Eq. (2.5). As our gluonic equation did not respect multiplicative 
renormalisability, we were not able to do so, and got an optimum as ~ 1.4. To see what 
the effect of the coupling on the quark propagator is, we shall vary as independently of Z^ 
although the two are really correlated. 

We show the result of this exercise in Eig. 4. As the value of the coupling grows, 
the quark propagator is enhanced near the origin. As as grows further, the propagator 
experiences oscillations, until the chiral solution is lost. It is interesting to note that the 
value of as that we obtained for the gluon propagator is very close to the critical value 
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Figure 5: Same as Fig. 4, but for with a transverse vertex multiplied by 






beyond which the chiral solution disappears. We also see that the criterion of Section 5.2 
can in practice work backwards, i.e. it not only predicts when there will be chiral solutions, 
but also when these will disappear. 

In order to investigate the sensitivity of the solution to our choice of the function T> 
in Eq. (3.11), we consider two possible choices for that function. In Fig. 5, we show the 
behaviour of our solution for for T^ = q"^ -\- p'^ , instead oi T> = [q^ + p^Y (as in Fig. 3). As 
expected, the two solutions have the same high-(j'^ behaviour and are equal at q^ = fi^. We 
see that the low-momentum behaviours vary by a modest amount, and that the solutions 
only differ in the region q^ sa //j, and only when the integrand oscillates a lot. This is 
where we do not trust our approximation for the angular integrals and hence the result is 
stable where our approximations hold. 
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Figure 6: The dashed curve shows the value of Uc for various functions r] (see text), and the plain 
curve shows the value of the intercept of the propagator for as = 0.75 Uc- 



Finally, Fig. 6 illustrates the insensitivity of our results to the choice of transverse vertex. 
We consider the functions rj = [[q"^ + p'^) j {q^ + p^Y] and plot the results in terms of n. We 
show the variation of the critical value of as, which changes by a factor 2 when n changes 
by a factor 100. We also show the value of the intercept at the origin of the propagator 
calculated at as = 0.75ac, and which hardly changes with n. This clearly illustrates that 
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the results we have obtained hold for a wide class of transverse vertices. 



6. Conclusion 

We have shown that the transverse part of the vertex plays an essential role in the quark 
DS equation in the axial gauge. It restores multiplicative renormalisability, allows chiral 
symmetry breaking and provides solutions which match with perturbation theory at large 
q^ . Hence we have obtained an equation that possesses all the properties that are required 
in QCD. 

The solution that we obtain in the chiral phase has a pole at (j'^ = despite the fact 
that we expect quarks to be confined. This contrasts with the results of a similar analysis 
carried out in ref. [15] for the gluon propagator of ref. [8] in which it was shown that for the 
confining propagator of ref. [8] the quark propagator does not have a singular behaviour at 
(j'^ = and can therefore be considered to be confined. The present formalism can be used 
to solve for the chirality breaking solution, as we shall explain in a future publication [16]. 
It is an interesting question to see whether these solutions are confined, and hence whether 
confinement and chiral symmetry breaking are related. 

It has not so far been possible to obtain a solution to the (modified) DS equations for 
the gluon propagator for time-like momenta. This is because we expect that in this regime 
the propagator develops an imaginary part with cuts corresponding to the thresholds for 
glueball production etc. and consequently such a solution involves coupled non-linear equa- 
tions for the real and imaginary parts of Zji^k^). Therefore the corresponding solution for 
the quark propagator we have obtained is only valid for quarks propagating with space-like 
momenta (in the t-channel). It is tempting to try to analytically continue the solution ob- 
tained into time-like momenta by fitting the solution obtained to known analytic functions 
and possibly then using the DS equation for the quark propagator to extract information 
about the gluon propagator with time-like momenta. Unfortunately this has not been pos- 
sible. We have obtained several very accurate fits (errors nowhere worse than 2% ) to the 
curve shown in Fig. 4, using different parametrisations inside various analytic functions. 
Whereas these all fit the quark propagator remarkably well in the fit region they give wildly 
different projected behaviours for the quark propagator for time-like momenta. We there- 
fore have to accept for the moment that the important problem of describing the exchange 
of soft gluons and quarks in the s-channel remains unsolved. 

As we explained earlier, the equation for the gluon propagator also will have to be 
modified to incorporate a transverse vertex and recover multiplicative renormalisability. 
Furthermore, as has been pointed out in ref. [19] the DS equations for the quark and gluon 
propagators are coupled and after having obtained the solution for the quark propagator 
one must check the assumption that one may neglect quark loops in the equation for the 
gluon propagator. In ref. [19] it was shown that for the gluon propagator of ref. [8] this was 
not the case and the gluon propagator had to be modified accordingly. A similar analysis 
should also be carried out for the gluon and quark propagators discussed here. 
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7. Appendix: Kernels 

The kernels entering Eq. (4.1) are given by the following expressions, where we define 
p = P/q"^. 

For the BZ equation (3.6): 



A.(.) = ^''(.-4-i(^J 



(7.1) 



and 



AM = <,(l-,)^ + .(,-l)(i(l + i)-^ 



+0(.-4)(l+,)^l-i(f+l 



(7.2) 



For Fq. (3.15), and for rj a function of 



Jt£l 



1 + P 
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(7.3) 
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[/>3 + p^ + 18/> + 16] 
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2p 
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